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Abstract 

We discuss the difference between the distribution of secondaries measured in terms 



Oh 

, of pseudorapidity and that using the correct rapidity variable. We show a set of examples 

obtained using Monte Carlo simulations. We also consider the production of particles of 
low transverse momentum where coherence effects may occur, which are not yet included 
in the present Monte Carlos. 
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The purpose of this note is to recall the difference between rapidity, y, and pseudorapidity, 
77. This is not a new topic pQ. However, since the results of the LHC detectors are shown 
mainly in terms of pseudorapidity, it is timely to discuss effects which may be caused by the 
^ ! identification of rj with y. 

1 Definitions of rapidity and pseudorapidity 

The Lorentz-invariant single-particle inclusive cross section has the form 



d 3 a „ d 3 
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/(AB^CX) = E— = E (1) 
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where E, pi and pt are the energy, longitudinal and transverse momentum of the outgoing 
particle C with respect to the incoming proton direction; and d 3 p/2E is the Lorentz-invariant 
phase space 

H = d*p5(p 2 -m 2 )Q(E), (2) 



/(AB^CX) = -_=, (4) 



where m is the mass of particle C. For high-energy hadronic collisions, experiment shows that 
the longitudinal momentum may take almost any value, whereas the transverse momentum is 
usually small. 

The rapidity y is a convenient way to display the pi dependence of the data, where 

1 ( E + p L \ ( E + Pl \ 

y = -log = log , 3 
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where = m 2 + p\. The value of y depends on the reference frame, and, at first sight, it 
looks like a rather ungainly construct. However, there are advantages to present high energy 
data as a function of rapidity. First, we see dy = dp^/E, and so the Lorentz-invariant cross 
section 

d 3 a 
ndydp\ 

Second, rapidity is additive under Lorentz boosts along the beam direction. For example, under 
a boost of velocity u 

E — >■ l(E + upi) , p L -+7(p L + uE), (5) 
where 7 = 1/ Vl — u 2 . Hence 

!, f 1 + u \ 

y -»■ y + 2 g \r^lJ = ^ + ^ boost - ( 6 ) 

Thus the difference in rapidity of two particles is not changed by boosts along the beam axis. 
In the non-relativistic limit, when the boost t; < 1, we have E — > m, p — > mv and thus the 
rapidity y — > velocity v (hence its name). 

For a massless particle E — \p\. Then the rapidity can be directly expressed in terms of the 
scattering angle 9 

In general, we define pseudorapidity as 



1 /l + cos0\ , 9 

V = 2 lo HT3-^J = "log tan-. (7) 



Q 

r] = -log tan-. (8) 

Clearly, for massless particles the rapidity and pseudorapidity are the same. However, this is 
also the case for particles with pt 3> m and >m. 



2 Expectations at the LHC 

We now come to the single-particle rapidity distributions measured at the LHC. From the 
theoretical viewpoint, both the phase space of particle C and the matrix element should be 
written in terms of covariant variables, that is in terms of y and pt- Long ago, Regge theory 
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Figure 1: A sketch showing the inclusive production of particle C. 



predicted a flat rapidity distribution in the central region. Indeed, in Regge theory it assumed 
that pt is limited and (pr) does not increase with energy. At very high energies, the inclusive 
cross section is described by the emission of secondaries from the Pomeron. The probability 
of the emission of particle C is given by the effective Pomeron-C vertex, (3pc, which does not 
depend on rapidity (2J [3] , see Fig. CD We have 

d& n2 o , \ a(0)-l a(0)-l / n x 

If y is measured in the pp centre-of-mass frame, then 

Sacbc = m p m T exp(±y + Y/2) (10) 

(where sac = (pa + Pc) 2 , Y = \n.{s/mi) and a(0) is the intercept of the Pomeron trajectory), 
and the independence of the distribution on y is explicit. The cross section may also be written 
in terms of parton distributions 

da f da(gg C) ( 1 x x 



dydp\ 



with x 12 = (m T /y/s)exp(±y). If the gluon distributions can be described by a fixed power, 
g(x) ~ x~ x , then again we see that the single-particle multiplicity distribution is expected to 
be flat in y. Indeed, the present Monte Carlos are based on ( ITTi) . and yield a flat distribution 
in the central region, see Fig. |2j The PYTHIA 6.4 Monte Carlo (ATLAS) |4J was used to 
demonstrate the effects in Figs. EH3 



3 The effects of approximating y by r] 

The single-particle distributions so far measured at the LHC have limited particle identification. 
The three-momentum of the particle is measured, but not, in general, its mass. Without 
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Figure 2: The pr- rapidity correlation obtained for pions from a Monte Carlo sample. 



knowledge of m it is not possible to calculate y. Therefore the results are presented in terms of 
the pseudorapidity rj. For 77 = y we need px 3> m. However, for low pt the equality does not 
hold. In particular, a particle with low pt may have small y, but rather large 77, since for low 
Pt this slow particle tends to be emitted at a small polar angle 9. For simplicity, we consider 
Pl ^> m. Then the difference 



The log may be quite large as pr — > 0. The most interestingregion, 77 ~ 0, in the multiplicity 
distributions measured at the collider, contains a singularit}o at 77 = 0. As a consequence the 
low pt particles are spread over a larger 77 region, depopulating the 77 ~ region. This is clearly 
seen in Fig. [3] as a dip in the 77 distribution for 77 ~ 0. For distributions corresponding to larger 
values of pt the dip is less pronounced. 

This trivial kinematic effect produces some 77 dependence of (pt) measured in the central 
detector. Since low pt particles are spread over a larger number of low 77 bins, the measured 
value of (pt) decreases with increasing \rj\, see Fig. HI This effect is not present in the y 
distributions. 

Of course, the majority of secondaries are pions. Therefore the multiplicity distributions 
may be presented in terms of y by assigning all particles to be pions. Nevertheless there are 
two sources of potential difficulty First look at Fig. [5j It shows the fractions of different sorts 

^or fixed pt, the Jacobian between the dy and drj bins is drj = (E/\p\) dy. For low pt, it has a singularity 
at \p\ -> 0. 




(12) 




Figure 3: The pp-pseudorapidity correlation obtained from the same Monte Carlo sample as Fig 12 
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Figure 4: The mean value of p T as a function of pseudorapidity (upper points) and rapidity (lower 
points). 
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Transverse momenta, p T , MeV 

Figure 5: The Monte Carlo prediction for different types of minimum bias particles as a function of 
their p T . Although the pions dominate, there are a appreciable number of soft electrons in the very 
lowest px bins and a significant kaon and baryon component at the larger p T values. 

of secondaries according to a PYTHIA Monte Carlo simulation. First, we see the presence of 
electrons (arising from Dalitz pairs) in the very lowest px bins. If the electron is assigned the 
mass of the pion, then its 'effective' rapidity becomes close to zero, which produces a prominant 
erroneous peak for y ~ 0. Second, at larger p?, the kaons and baryons, which are misassigned 
as pions, lead to an erroneous dip for y ~ 0. This dip disappears as we go to larger values 
of px- Note however, that to neglect the difference between y and rj, we need a value of pr 
larger than the mass of the corresponding hadron. All these effects are seen in Figs. E] and [71 
Note, from Fig. [51 that in the larger pt region, the fraction of kaons and baryons is appreciable. 
Only about 60% of secondaries are pions. Recall also that the region of low px and small r\ is 
very poorly described by the present Monte Carlo models. We should therefore take care when 
discussing this domain. 

4 Examples of the necessity to use y rather than r] 

To demonstrate the importance of using the correct rapidity variable, we consider the coherence 
effects in low pr particle production. In the Regge approach, the p T distribution of secondaries 
produced in the central region (y ~ 0) does not depend on energy, but reflects only the internal 
Pr-structure of the Pomeron. Therefore we expect that the particle density dN/ dydp T measured 
at low pt will increase with energy just like the density dN/dy integrated ovecl over p?. That 

2 The region of a large px gives a small contribution. 
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Figure 6: The pr-rapidity correlation assuming that all the particles of Fig. [5] have the mass of the 
pion. 




Figure 7: The same as Fig. El but omitting the pion component and taking the other particles to 
have the mass of the pion. 
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is, by about 60% in going from the collider energy of 0.9 to 7 TeV. In terms of Monte Carlo 
models such an effect is simulated by the number of Multiple Interactions (MI), which increases 
with initial energy of the colliding protons. 

On the other hand, such an approach does not account for the coherence effects which 
may be quite important for low momenta secondaries. For low pt particles, interference arises 
when a particle is radiated from different lines of Feynman diagrams which make up the whole 
matrix element. One example of this is soft gluon emission. Another example is the Hanbury- 
Brown— Twiss effect, where two identical particles can be radiated from different places in the 
interaction region. Their interference leads to a peak in the distribution of the two identical 
particles when their momenta are close to each other. The width of the peak can be used to 
estimate the size of the interaction region [5]. 

First we discuss soft gluon emission, which leads to the so-called limiting energy behaviour. 
The radiation of soft gluons from several quarks and/or gluons develops in a coherent way. The 
reason is that different sources of secondaries should act as a single source with an effective 
colour 'charge' equal to the 'vector^ sum of all the colour 'charg gluon of large wavelength 

cannot resolve these smaller details. Recall that, due to the Lorentz time dilation, the low 
momentum hadrons are formed at the initial stage, before the fast particles [6]. These fast 
particles, which are formed later on, fly out of the interaction region and leave the slow hadrons 
intactEl- Such a limiting behaviour of soft secondaries was first predicted, and observed, in e + e~ 
annihilation E], and then extended to the pp case [9J. 

Also, note that, if the incoming hadron represents a quark-antiquark (meson) or quark- 
diquark (nucleon) system, then the t-channel multi-gluon exchange corresponds dominantely 
to a colour octet. The same colour-octet exchange dominates the BFKL amplitude. So even in 
the case of a large number of gluon exchanges (or a large number of MI), the low pt secondaries 
will be produced mainly via radiation caused by colour-octet flow. 

Thus, in the ideal case, we may expect a limiting energy behaviour of the soft particle 
density. That is, at very low p?, the value of dN j ' dydp\ should not depend (or depend only 
weakly) on the initial energy, see Ref. jH]. However, as was shown above, in order to measure 
the density of low-pr particles at small y, it is important to use the correct 'rapidity' variable, 
and not just the angle (that is, rf). Otherwise the expected 'limiting' behaviour may be affected 
by an admixture of other sorts of particles, some lighter and some heavier than the pion. 

The Hanbury-Brown— Twiss correlation can be viewed as the first step in the formation of 
condensates of bosons which arise from Bose-Einstein statistics. At high energies when the par- 
ticle density increases, and the number of particles exceeds one in the elementary Ax Ap cells of 
some domain of configuration space, there may be two sorts of such collective phenomena. One 
is of pure classical, thermodynamic origin, like elliptic flow in quark-gluon plasma. Another 
one is based on quantum mechanical interference. We consider the formation of some mul- 
tiparticle coherent state (such as in superconductivity or in a disorientated chiral condensate 

3 An analogy is a supernova explosion where the remnants of the star are left intact. 
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[10J). The coherent emission of soft gluons by the vector sum of the whole colour charge may 
be considered as the first signature of the multiparticle coherent state. Since the wavelength of 
low pt particles is large, it is expected that such collective phenomena will occur, and be best 
observed, in the low pt region. In comparison with lower energies, these collective, coherent 
effects should be enhanced at the LHC due to the larger density of low pt secondaries. 

In summary, to observe such phenomena we need a careful study of very soft particle 
distributions using the correct Lorentz covariant variables. 
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